A comparative study of crumpling and folding of thin sheets 
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Crumpling and folding of paper are at first sight very different ways of confining thin sheets 
in a small volume : the former one is random and stochastic whereas the latest one is regular and 
deterministic. Nevertheless, certain similarities exist. Crumpling is surprisingly inefficient: a typical 
crumpled paper ball in a waste-bin consists of as much as 80% air. Similarly, if one folds a sheet of 
paper repeatedly in two, the necessary force becomes so large that it is impossible to fold it more 
than 6 or 7 times. Here we show that the stiffness that builds up in the two processes is of the same 
nature, and therefore simple folding models allow to capture also the main features of crumpling. 
An original geometrical approach shows that crumpling is hierarchical, just as the repeated folding. 
For both processes the number of layers increases with the degree of compaction. We find that for 
both processes the crumpling force increases as a power law with the number of folded layers, and 
that the dimensionality of the compaction process (crumpling or folding) controls the exponent of 
the scaling law between the force and the compaction ratio. 

PACS numbers: 46.25.-y, 46.32. +x, 46.70.-p, 62.20.-x 



It is easy to verify that the maximum number of times 
one can fold a sheet of paper is only 6 or 7, which surpris- 
ingly is independent of the initial size of the sheet. Quan- 
titatively, elasticity theory allows to write the relation 
between the compaction force and the number of times 
one can repeatedly fold a piece of paper in two. This fol- 
lows from the scaling of the bending rigidity B with the 
thickness h of the folded sheet B = Eh^ / 12(1 - u^), 
where E is the Young modulus and u the Poisson ratio. 
For a sheet of initial size Di x D2 folded along the direc- 
tion the energy, Ec^ injected in the system should be 
compared to the typical energy dissipated in the fold. Ec 
can be readily written as Ec = FD2, where F is a char- 
acteristic compression force applied along the direction 
Di. Since most of the folded sheet remains flat and the 
region which is irreversibly deformed is straight along D2 
(i.e. its gaussian curvature equals to zero ; See Fig. [T]), 
the energy dissipated in the fold can be estimated from 
the elastic bending energy, E^i, concentrated in a region 
of length D2 and width h with a curvature [1] ; this 
leads to Eei ~ BhD2/h^ = BD2/h. The balance of these 
two energies leads to F ^ B/h. When the sheet is folded 
n times repeatedly leading to the hierarchical creation of 
folds, its effective thickness and bending rigidity become 
hn 2^h and Bn 2^^B assuming no slip between lay- 
ers (this hypothesis becomes increasingly consistent for 
large folding events n). Consequently, for a sheet folded 
n times, the energy balance gives 

F{n) ^ Bn/K - Fo2^^ , (1) 

where Fq ^ B/his the elementary force needed to create 
a unique fold. Thus the force is independent of the initial 



size of the sheet and grows exponentially with the number 
of folding events n. The exponential dependence is the 
reason why one cannot fold a sheet indefinitely by hand or 
by applying a finite force; a simple estimate can be made 
taking typical values E = lO^Pa, h = 10~^m, leading 
to a prefactor Fq « 1 N. Then, F becomes of the order 
of kiloNewtons already for n = 6, which is larger than 
the maximal force any person can exert and then sets a 
limit on the number of successive foldings that can be 
achieved. 
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FIG. 1. Hierarchical folding of a sheet in different dimension- 
alities, (a) a Id-like sheet folded in Id; (b) a Id-like sheet 
folded in 2d; (c) a 2d- like sheet folded in 3d. The three types 
of folding processes are referred to as Id, 2d and 3d com- 
pactions. 

Repeated folding in two is of course not the only pos- 
sible way to fold; here we consider three basic regular 
processes (Fig.jl]), which can be seen as prototypical fold- 
ings in various dimensionalities. In cases (a) and (b), the 
sheet can be thought of as a Id-like sheet, since it is 
folded along one direction only. However compaction of 
case (a) is not isotropic, contrary to case (b). The lat- 
ter can be seen as isotropic compaction within a disk, 
whereas the former can be seen as unidirectional com- 
paction inside an elongated rectangle such as in [2 . Fi- 
nally, for case (c) the sheet is a truly two-dimensional 
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object that is compacted in a sphere. The number of 
folded layers N = hn/h after n folding events and the 
related compaction ratio (j) (defined as (j) = I^/A, where 
D and A are the initial and final size of the sheet) then 
depends on the precise geometry and dimensionality of 
the compaction process. Folding a (i- dimensional sheet 
in {d + 1) -dimensional space, N = 2^, whereas for the 
Id-ld case (a), = n. Also, for foldings (a) and (b), 
N = (j)^ whereas in the Sd-like folding (c). A/" = (/)^. Using 
similar arguments leading to Eq. ([T]), one finds a generic 
power law relation between the force F, the compaction 
ratio (j) and the number of folded layers N: 
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Id compaction, case (a) 
2d compaction, case (b) 
3d compaction, case (c) 



F{N) = FoN - Focj) (2) 
F{N) = FoN^ - Fo(^^ (3) 
F{N) = FoN^ - Fo(l)^ (4) 



It is important to note that we described folds as the re- 
sult of an irreversible process occurring in a small region 
of size h and characterized by a zero gaussian curvature. 
Consequently, energy scalings are different from those ob- 
tained for singular ridges p| ^ for which irrreversible 
process are not taken into account. Moreover, our esti- 
mate for the plastic energy dissipated in the fold should 
be taken as a lower bound because the contribution of 
vertices [5] and other possible length scales associated 
with plastic events have been neglected. Our approach 
is inspired by crack modelling in the framework of linear 
elastic fracture mechanics where the process zone near 
the crack tip is neglected and the dissipation is described 
by a single parameter : the fracture energy. 

One may wonder to which extent our models for reg- 
ular folding describe our crumpling measurements and 
whether crumpling process can be viewed as arising from 
successive folding events. To compare the regular folding 
with the crumpling of paper, we first show experimen- 
tally that here the force also increases as a power law 
with the degree of compaction, with exponents in ac- 
cordance with predictions dependant on the dimension- 
alities of the compaction process. For this purpose, a 
sheet of paper of characteristic size is placed into a 
rigid cylindrical cell, in which a piston connected to a 
force transducer that compacts the sheet into a pancake 
at constant speed (Fig. |2|i). The materials used here 
(Kraft and regular printing papers) have been chosen 
for their low ductility. The experimental force-distance 
curves show a very strong increase of the force upon com- 
paction (Fig. |2]3), for both 'virgin' sheets (sheets crum- 
pled for the first time) and 'trained' sheets (sheets crum- 
pled for the second or third time). For the latter the 
force-compaction ratio curves turn out to be independent 
of the initial preparation of the sheet inside the cylinder 
within the experimental accuracy. Then, the measured 
curves for different types of paper of similar properties 
(thickness h = 10 jam and Young modulus E lO^Pa), 
different sheet sizes, cells and compaction speeds can all 
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FIG. 2. (a) Schematic of the setup used to measure the 
force F{A) during crumpling compaction, (b) Typical force- 
distance curves obtained from crumpling experiments, with 
Kraft paper of size D = 15cm in a cylindrical cell of diame- 
ter 6cm. The three curves result from the 1** (blue circles), 
2'^^ (green squares) and 3^"^ (red diamonds) crumpling rounds 
(only 1 out of each 100 experimental points are drawn for 
clarity). Inset: Same data as a function of compaction ratio 
F(A/D) shifted with respect to each other for clarity on a 
log- log scale. The lines are the fits to a power law (Eq. [5]): 
the fitted values (a, /3) are (3.74 N, 1.1), (1.90 N, 1.4) and 
(1.43 N, 1.5) for the 1"*, 2^^ and 3^^ crumpling rounds re- 
spectively, (c) Exponent /3 of the power law fit (Eq. [5]) as 
a function of the paper size D for the 2^^ (green) and 3^^ 
(red) crumpling rounds. Inset: The probability distribution 
function of /3. (d) Characteristic force scale a of the power 
law crumpling force as a function of the paper size D. Inset: 
The probability distribution function of a. For (c) and (d) we 
averaged over many experimental realizations (Kraft and reg- 
ular printing paper of similar properties) for constant values 
of D; the error bars give the spread. 



be described by a power law: 



F(A) = a 



D 



(5) 



where a is a characteristic force scale and A is the gap 
between the piston and the bottom of the cell (Fig. |2|i). 
While the data range for A is small, this behaviour is ro- 
bust over all the 150 realizations. The statistical test 
for goodness-of-fit confirms the relevance of the power 
law in comparison with other fits. The exponent p as- 
sociated with the power law divergence is found to be 
/3 ~ 1.3 (Fig. |2]3), a value between 1 and 2, those ex- 
pected for ordered folding in Id and 2d. We argue below 
that this is due to the anisotropy of the compaction pro- 
cess in our experiment. Effectively, compaction here is 
quasi Id, since loading is applied mainly in one direc- 
tion. However, the setup also allows for compaction in 
the perpendicular direction, which would rather be a 2d 
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process. Moreover, we find the characteristic force scale 
a to be independent of size D and equal on average to 
2 N (Fig. [2]i) which is of the same order of magnitude as 
the characteristic force, Fq ^ IN, calculated above for 
the folding. 




FIG. 3. (a) Picture of a crumpled cross-section; the white 
lines represent two orthogonal directions used to extract the 
number of folded layers N . (b) Picture of a piece of an un- 
folded cross-section and its segmented edge. Scale bars are 
10mm. (c) Number of folded layers A/" as a function of com- 
paction ratio D/A for crumpled balls of Kraft paper (circles) 
and regular printing paper (triangles). The solid line is the 
curve N — {D//S)'^. (d) Probability density function p{l) of 
the lengths of segments on a log-log scale for a Kraft paper 
sheet crumpled in a ball of diameter A ^ 80mm and com- 
paction ratio 6. p{l) is compared to a log-normal distri- 
bution (continuous line) and a Gamma distribution (dashed 
line) with the same mean and variance as the empirical data. 
The error bars 51 and 5p oi the experimental pdf's are given 
by the bin width 31 and the estimated standard deviation 
8p — pI \fn of the corresponding histograms n(^). 

A second step towards understanding the analogy be- 
tween crumpling and folding is to establish the relation 
between the degree of compaction and the number of 
folds for the experimental crumpled configurations. To 
achieve this, we characterize the geometry of the crum- 
pled paper, using an original approach that makes use 
of the properties of folds and facets in cross-sections of 
crumpled samples. Sheets of different paper types and 
sizes D are crumpled into hand-made balls at different 
degrees of compaction. A cross-section is obtained by 
cutting the crumpled ball in two with a slowly moving 
hot wire [6]. The overall size of the crumpled configura- 
tion, A, is defined as the largest diameter of the result- 
ing cross-sectional area. In this cross-section, the number 
of paper layers is measured in two orthogonal directions 
passing through the center (Fig. |3|i), and subsequently 
averaged to obtain the mean number of folded layers N 



in the crumpled configuration. By this method, we en- 
sure that is defined as in the folding model introduced 
above {N = hn/h). The number of folded layers N can be 
described by a power two dependance N {D / /S.^ = cjp' 
on the degree of compaction, with a prefactor close to 1 
(Fig. [3]:). This important result is exactly the same as 
that observed for 3d folding (Fig. [l]^) showing that the 
geometry of folding and crumpling is the same. 

A possible difference between the two situations is that 
while the repeated folding is a hierarchical process, this is 
not at all clear for crumpling. To investigate whether the 
crumpling is also hierarchical, we characterize the lengths 
of folds and facets in cross-sections of crumpled samples. 
For this purpose, the cut crumpled sheet is reopened care- 
fully. Several uncrumpled pieces, with possibly several 
holes, are obtained and scanned (Fig. jsj)). The edges 
of their boundaries and holes are detected automatically 
and broken down into segments delimited by kinks [7], 
by using a 'Split and Merge' algorithm |8| for the seg- 
mentation. The planar two-dimensional cross-section of 
the crumpled sheet bears information on the full three- 
dimensional crumpled configuration: the ensemble of seg- 
ments samples the facets delimited by folds, so that its 
length distribution can be related to a characteristic dis- 
tance between folds or equivalently to the characteristic 
size of the facets. We can use this to asses the nature 
of the crumpling. To do so, we compare the distribu- 
tion of lengths p{t) with a log-normal distribution and 
a Gamma distribution; the former characterizes a hier- 
archical process [H [71 [9] , whereas the latter accounts for 
random processes [H [71 [9] . More precisely, a log- normal 
distribution describes a fragmentation process in which 
all pieces are broken successively into two parts, such that 
any new fragment is further broken into two pieces where 
the breaking point being uniformly distributed along the 
fragment [71[T0l|TT]. In contrast, the Gamma distribu- 
tion emerges from a fragmentation process where all the 
breaking points are uniformly distributed along the un- 
broken line, prior to the breaking that happens simul- 
taneously for all points. Fig. [3]i shows that both dis- 
tributions reasonably well describe the rapid decay of 
the tail of the distribution, but the Gamma distribution 
seriously overestimates the probability density at small 
lengths. This originates from the fact that a hierarchi- 
cal fragmentation process tends to generate less small 
fragments than a random one. A more rigorous test is 
done through the statistical ^^st for goodness-of-fit, 
which confirms that the log-normal describes better the 
data. We checked that this description is robust with 
respect to the chosen value of the threshold used in the 
segmentation procedure. The log-normal distribution ac- 
curately describes all the experimental data sets, and so 
the crumpling is hierarchical rather than random. Ear- 
lier simulations [12] of crumpled sheets and experiments 
on unfolded sheets [13] found a similar agreement with a 
log-normal distribution. 
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TABLE I. A summary of various results for the power- law 
variation of the force F oc (D/A)^ for different cases of the 
compaction of an x-dimensional object in (x + l)-dimensional 
space. We provide measured values of /3 from the literature 
and our theoretical prediction /S'^ of the crumpliness expo- 
nent. The latter follows from the dimension that is given 
by the geometry of compaction. For some cases, topologi- 
cal constraints and the material properties also influence the 
value of the exponent. The first set corresponds to results 
that can be approached by a mixture of case (a) and (b) of 
Fig. [l] The second and third sets correspond to cases (b) and 
(c) respectively. 



Crumpled Object and Ref. 


X 




/3 




Paper (this work) 


2 


1 < x"^ < 2 


L3 


1< < 2 


Mylar TO 


2 


1< X* < 2 


L89 


1< < 2 


Tethered Membrane [l7] 


2 




L85 


/^"^ ^ 2 


Rods [18] 


1 


1 


2 


2 


Rods 


1 


1 


2.05 


2 


Linearly Elastic Sheet [l2j 


2 


2 


4 


4 


Aluminum Foil [20^ 


2 


2 


5.13 


6 


Phantom Sheet [l2] 


2 


2 


2.66 


2.5 



The conclusion is that folding and crumpling are very 
similar in nature and the crumpling process can be 
viewed as arising from successive folding events. For 
ordered folding, simple models allow for predictions of 
the relations between force, compaction ratio and num- 
ber of folds. Surprisingly, these are found to capture 
the main properties of crumpling also, in particular the 
hierarchical structure of the folds and the power law re- 
lation between the force and the compaction ratio. The 
analogy with folding then allows to define the 'crumpli- 
ness exponent' P for various forms of crumpling process. 
Previous experiments and simulations in the literature 
have reported such exponents (3 for the power-law de- 
pendence of the force on the compaction ratio. It follows 
that also these cases can be explained in a satisfactory 
way using our arguments, i.e., solely by considering the 
dimensionality of the compaction process, the topological 
constraints and the mechanical properties of the material 
(e.g. ductility). Table [l| summarizes various exponents /3 
found in the literature, detailed below. 

For the first set of data, Matan et al. [16] used a com- 
paction set-up similar to the one used here, and found an 
exponent of /3 ~ 1.89 (this is the inverse of their exponent 
a). The aspect ratio of their cylinder (height/diameter) 
is much smaller than ours; we thus anticipate that the 
compaction is more 2d in nature, and hence one would 
expect a crumpliness exponent closer to 2 than in our ex- 
periment, which is indeed observed. The value of this ex- 
ponent can again be understood as a compaction process 
lying between Id and 2d. As our arguments are based 
on dimensionality, they allow to predict only bounds for 
this type of experiments. The simulations of [17 on teth- 
ered membranes that are compacted found a value of 
P ~ 1.85. Except that loading is now biaxial, the com- 
paction process is in fact similiar to that of case (a) since 



the "height" fluctuations of the membrane are small. If 
the two directions were independent we recover a number 
of folds given by ex and thus force F = NFq cx 
However, folding in one direction is inhibited by folding 
in the other direction. This effect will decrease the to- 
tal number of folds leading to a crumpliness exponent 
/^"^ ^2. 

For the second set, the analogy with case (b) is com- 
plete and both experimental [19l, theoretical [18 and nu- 
merical [19 results are in agreement with the simple ar- 
gument for hierarchical folding predicting a crumpliness 
exponent = 2. 

The third set deals with experiments and simulations 
of 2d sheets crumpled inside 3d spheres. The linearly 
elastic sheet [12] is a perfect example of case (c) for which 
the crumpliness exponent P'^ =4, in agreement with the 
simulations. The aluminum foil represents ductile sheets 
with plastic deformations [20] and the phantom sheet is 
a sheet that can cross itself in the simulations [12 . These 
are somewhat more complicated cases; however we can 
provide an estimate of the exponent P for these two cases 
also. For the aluminum foil, because of their ductility 
one has to modify the estimate of the elastic energy E^i 
of a folded sheet. Again, most of the folded sheet re- 
mains flat but the elastic energy is now concentrated in 
a region of length D and width l//^c with a curvature 
f<ic^ which is a material constant: the curvature scale at 
which the material yields. This is the main difference 
with the purely elastic case for which the thickness h 
sets the scale over which the sheet is curved. Using the 
same arguments, the balance of the compaction energy 
and the elasto-plastic energy for the plastic case then 
leads to F ^ Bntvc ^ Eh^tZc As N = hn/h = (/), 
one consequently finds F{N) ^ N^Fq ~ i^o^^ leading 
to an estimated value /^"^ = 6, which is in fair agreement 
with experimental results [20] . For the simulations of the 
phantom sheets, the absence of steric interactions implies 
that F{N) ~ NFq which is similar to the Id case. The 
number of folded layers N is then related to the com- 
paction degree through N ^ D'^ / < S >^ ^l^f^ where 
< 6* > is the mean facet size and V/ ~< 6* > is the 
average volume occupied by the sheet. For high com- 
pactions, it is known that ~ F^g^ where Kg is the 
radius of gyration and (ij is the fractal dimension [21j . 
Thus we find F cx A/" - (\)^\ with p"" = df ^ 2.5 EHIH]. 

Finally, these arguments allow also to explain why a 
wastebucket fills up so quickly when waste paper is crum- 
pled into a ball. Using the equivalence between crum- 
pling and folding, the wasted volume AV/V can be esti- 
mated from the folded case. In the 3d case, this is given 

by 

A3 =^-^''^'- (6) 

For N = 2^ and typical paper {h = 10~^m, D = 0.2m), 
one has AV/V ~ 75%, which is an excellent estimate for 
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the experimental observation that crumphng is a very 
inefficient compaction process. In conclusion, the obser- 
vations presented here demonstrate a non-trivial relation 
between the force of compaction and the geometry of the 
crumpled configuration. A potential application of this 
result would be to invert this problem, and deduce the 
force through analyzing cross sections of crumpled sheets. 
Since the arguments presented are generic, they should 
hold also at the nano-scale and could provide a simple 
framework to understand crumpled graphene structures, 
such as graphene-based supercapacitors [Mj [15] . 
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